Field-induced Kosterlitz-Thouless transition in the N = Landau level of graphene 
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At the charge neutral point, graphene exhibits a very unusual high resistance metallic state 
and a transition to a complete insulating phase in a strong magnetic field. We propose that the 
current carriers in this state are the charged vortices of the XY valley-pseudospin order-parameter, 
a situation which is dual to a conventional thin superconducting film. We study energetics and the 
stability of this phase in the presence of disorder. 

PACS numbers: 72.10.-d,73.21.-b,73.50.Fq 

The initial experiments of the quantum Hall effect 
(QHE) in monolayer graphene discovered the quantum 



Hall plateaus a xy 
ing factors v 
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4 (N + \) (e 2 /h) (N e Z) at fill- 
±2, ±6, ±10,- • •§. Here, £ B 
is the magnetic length, p the carrier density measured 
from the charge neutral point, and the factor of 4 arises 
from the spin and valley (K and K') degeneracy. Re- 
cent experiments under stronger magnetic fields, on the 
other hand, showed additional plateaus at v = 0,±1,±4 
|^|, ||. Further experiments Q| at v — in high qual- 
ity samples revealed a rapid divergence of the longitu- 
dinal resistance R xx at a critical field B c . Interestingly 
such divergence fits the Kosterlitz-Thouless (KT)|J form 
Rxx ~ e a /VB a -B ovcr three decades of the resistance j|. 
Moreover, for B < B c the resistance saturates at low 
temperature to a value much larger than the quantum of 
resistance [|] , a behavior qualitatively different from con- 
ventional thermally activated transport in strong mag- 
netic fields. 

Since the critical field B c lies in the regime where the 
4-fold degeneracy of the N = Landau levels (LLs) is 
split, it is important to understand the cause of such 
splitting. There are several theoretical proposals. The 
common theme is the observation that the exchange ef- 
fect of the long-range part of the Coulomb interaction 
favors the spontaneous polarization of the real and/or 
valley-spins § | |, 0, 0, £| || gj, || g |§. 
(the valley-spin is a SU(2) pseudospin variable, which 
we denote by T; its z-component T z = +1 (—1) cor- 
responds to K (K'), respectively.) In one of the pro- 
posals, it is argued that the Zeeman energy favors the 
polarization of the real instead of the valley-spin jll], [l2|| . 
In another it is argued that the short-range part of the 
Coulomb interaction favors the spontaneous polarization 
of the valley-spin so that K and K' becomes unequally 
populated §, § § g § §. Because the TV = 
Landau states associated with K' and K localize on 
complementary sublattices (A and B, respectively) this 
amounts to a charge-density-wave (CDW) modulation 
which breaks the A-B sublattice symmetry. 

The type of diverging resistance observed in Ref. [0 is 



difficult to account for in the real spin polarization sce- 
nario since in this scenario there arc spin-filtered counter- 
propagating edge states that give rise to a metallic con- 
ductance of 2e 2 /h &tv = p] H, H|. Similarly in the 
CDW scenario ]13|, |b| [ll], |l6| it is difficult to explain the 
KT type resistance divergence, and the high-resistance 
metallic state below the critical field 0]. In addition to 
the above, there is a work (Ref. [pp| ) which claims an ex- 
plicit valley symmetry breaking term is consistent with 
the lattice point group symmetry. Like others, this work 
can not account for the KT behavior. 
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FIG. 1: (Color online) (a) The Kekule bond-density wave 
order with two defects marked by a filled circle. The defects 
are charged as they support a midgap electron state, (b) 
The U(l) phase (j> = tan -1 (T v /T x ) corresponding to the bond 
order pattern (a), (c) The valley-spin polarization ratio as a 
function of magnetic field B and sample mobility fi. (d) Same 
as (a) but for the pseudospin-supercurrent jj c /i^i cl ° an ' , where 
jz( clean) j g ^ e vame m -j-jjg c j ean ii m jt. The number of orbitals 
per valley and per spin is N$ = 50. 
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Motivated by the KT behavior and the highly resistive 
metallic state, wc propose an alternative scenario: the 
degeneracy splitting at v = is due to a spontaneous 
ordering of the pseudospin on the T x — T y plane (XY 
pseudospin ferromagnet). This involves a spontaneously 
generated hybridization between the N = LLs associ- 
ated with B (K) and A (iC), and is represented by the 
ground state wavefunction 
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is the creation operator for an electron in 
TO-th N = LL orbital at valley r = K, K 1 with real 
spin s =T)|- This type of order also breaks the lattice 
translation symmetry due to the mixing of B and A, and 
represents a bond-density-wave of some kind (Kekule or- 
der) [see Fig. 0(a)] @ [jj, II H Hi- The P hase <t> 
of this hybridization matrix element is the U(l) phase 
angle representing the direction in T x — T y plane, T = 
(cos^>, sin</>, 0), and associated with the sliding degrees 
of freedom of this density-wave. The low-energy charged 
excitations arc vortices and antivortices [Fig. 0(b)]. We 
study their binding-unbinding transition driven by mag- 
netic fields or disorder (see below). 

We now take a look at, piece by piece, the Hamiltonian 
for graphene in a strong magnetic field and the associ- 
ated energy scales, to address the plausibility of the inter- 
valley coherent state. (1) At the charge neutral point the 
LL separation is ^/2hv F /l B ^ 400a/bJtJ[K], which is the 
largest energy scale of the problem. Therefore in the rest 
of the paper we perform projection onto the N = LLs. 

(2) The Coulomb interaction He, which is the sec- 
ond largest energy scale of the problem, is approxi- 
mately symmetric under rotation in the combined space 
of real and valley-spins. The exchange energy is E^f ~ 
^/2(e 2 /d B ) ^ 120/B[T][K] §. 

As a result it favors the polarization of the SU(4) spin 
albeit it does not care whether the polarization should 
occur in the real spin, or valley-spin, or some combination 
of both§. 

(3) We now describe the parts of the Hamiltonian 
which break the SU(4) symmetry, 7?sb = / d 2 rHsB, 



-A Z S Z - U \S\ 2 - U Z T 2 - U X (T 2 + T 2 ).{2) 



Here S is the real spin operator. The first term in H$b 
represents the Zeeman energy where A 2 = gp^sB ~ 
1.3 x (B[T])[K]. The short-range part of the Coulomb 
interaction is not SU(4) symmetric and gives rise to Uq 
and U z . They can be estimated from the on-site and the 
nearest neighbor interactions, and are smaller than E(? 
by a factor a/^sH [|, where a is the lattice constant. 
While Uq favors the real spin polarization, U z favors the 
CDW phase (T z ^ 0)§ |, 0. 



On the other hand, U± term can arise from the 
electron-phonon interaction. One such example is the in- 
plane optical mode at K point, whose interaction with 
electrons can schematically be represented as[[l7| pi] 



H± = .f 



cfru-^ra^+^u 2 , (3) 



where u = (u x ,u y ) represents the (uniform) Kckulc- 
type distortion of the lattice, k measures the clastic 
energy, and Nc is the total number of carbon atoms. 
The two sets of Pauli matrices, {<J x ,y,z,o\ an d {T x ,y.z.o\, 
act on sublattice (A, B) and valley (K, K'), respectively; 
^{t x ± iT y )a x ip oc T x ± iT y oc e ±l * serves as the U(l) 
order-parameter of the Kekule bond density wave. Upon 
integrating out the phonon, this give rise to the U± term 
with U ± - 2.0 x (5[T])[K]. Note that U ± is comparable 
to Ajg. This can be traced back to the strong coupling 
between the K phonon and electrons p5|, |2^ |. Out-of- 
plane lattice distortion is studied in Rcf. Jl6|, and shown 
to contribute to U z which are much weaker than that as- 
sociated with in-plane modes in graphene [24|. The SU(4) 
breaking terms are summarized in Table 

Since the SU(4) symmetric part of the Coulomb in- 
teraction is much stronger than the symmetry breaking 
parts Hsb it is the former that sets the basic energy scale 
for the SU(4) symmetry breaking. The symmetry break- 
ing terms simply select the way the SU(4) symmetry is 
broken: they determine the nature of the ordered phase. 
Although Uo,U z ,U± in Table Q all have similar energy 
scales, it suggests the U(l) broken inter-valley coherent 
state (|l]) is a reasonable candidate for lifting the degen- 
eracy of the N — Q LLs. 

We now describe the field-induced transition at zero- 
temperature using the self-consistent Hartree-Fock (HF) 
theory. To account for competition between interaction 
and disorder effects, we allow the XY pseudospin order- 
parameter to be spatially inhomogencous. The matrix 
element of the HF Hamiltonian, in the Landau gauge, 
can be written in the form[p7|, 

(ma\H HF \m'a') = ^ Xm '^ v i%,x m -x m , 
i 

x (q) + U! a , (q) + (q) + <W U imp (g)] . (4) 
Here the system size is L x L — 2tt£ 2 b N^ 1 X m — 2ir£ 2 3 m/L 





residual symmetry 


energy scale 




[(spin) x (valley)] 






No x SU(2) 


1.3[K] x B[T] | 




Uo 


SU(2) x SU(2) 


1.0[K] x B[T] \ 




u z 


SU(2) x Z 2 (CDW) 


0.5[K] x B[T] J 






SU(2) x U(l) (Kekule) 


2.0[K] x B[T] |J 



TABLE I: SU(4) symmetry breaking terms, with the pattern 
of symmetry breaking and the energy scales. 
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(to = 1, 2, • ■ • , Nj,), and a = 1, . . . , 4 is the index for 
(real)spin and valley. The Hartrec and Fock potentials, 
and the anisotropic interaction are given as follows: 



2nt% 



V c {q)e- q -^Y.^"^ 



(5) 



L 2 

-Jh 
2tt 



cto 



3 ' 

is 



where Vc(q) = 2ne 2 /eq and T J 
The SU(4) order-parameter in Eqs. (|[ ||) 
determined sclf-consistcntly from A CT ' CT (qr) = 
V Ew(4ay CT ')e- i&Xm ^^ B> x m ,-x ? - The 
disorder potential £/i m p(<z) is given in terms of 
charged impurities located randomly at R/ by 



U imp (q) = L- 2 J2?L7e lqRl Vc(q)e 
The disorder strength is described by the impurity 



-lq 2 /i-iq x q v /2]ei 



imp 



(h/eB)n- im p, where n 



imp 



N imp /L 2 is 



hlling v\ 

the impurity density. In terms of the zero-field mobility 
of graphene, /i = a xx /ne ~ 2Oe//in; mp 0. the disorder 
strength is determined by a product B x /i. 

In the clean limit, Ni mp — 0, the ground state is fully 
pscudospin polarized due to the Coulomb exchange. In 
a weak field or in a dirty sample, on the other hand, 
the Coulomb interaction plays a minor role and the four 
(nearly) degenerated LLs are equally occupied. As shown 
in Fig. [I] (c), the pseudospin polarization ratio diminishes 
in the weak field and low mobility limit. In particular 
when B < B* ~ 10/(/4l0 4 cm 2 /Vs])[T] the pscudospin 
symmetry is restored by disorder. Here B* , referred to 
the dashed line in Fig. |l| (c), plays the role of the mean- 
field critical field. 

Even B > B* , vortices and anti- vortices tend to de- 
stroy the stiffness of the XY pseudospin order. An inkling 
of this KT transition could be seen in the unrestricted HF 
calculation as follows. In the pseudospin XY (quasi- long- 
range) ordered phase, angle-twisted states which have 
(j) = QX m in Eq. ([l]), namely pseudospin-supercurrent- 
flowing states, arc mctastable, because of the finite 
stiffness. Such states can be selected by artificially 
starting with an initial order-parameter Axu'iq) = 
<5q,Qx hi the initial step of the self-consistency loop. 
One can then monitor whether or not the pscudospin- 
supcrcurrcnt generated can survive as one iterates the HF 
calculation^. Charged impurities generate quenched 
vortex-antivortex pairs that randomize 0, and hence 
the pseudospin-supercurrent is expected to vanish in a 
weak field and low ps (see Fig. 0.) In Fig. |l| (d), the 
pseudospin-supercurrent of the metastable state, given 
by H| j Z sc = TjE q ,ss'r^q x i}V c (q)prs(-q)p-rAq), 
is plotted, where r = ±1 for K and K' , and p T s(q) is 
the density operator for valley r and spin s. Here note 
that pscudospin-supercurrents are charge neutral objects 



and thus are not related to true charge-currents. We 
calculate j* c to discriminate whether vortices and an- 
tivorticcs arc bound or unbound. As Fig. [I] (d) indi- 
cates, the pseudospin supercurrent drops around B c ~ 
40/(m[10 4 cm 2 /Vs]) [T] > B* . 

When B ^> B c and when the temperatures is suffi- 
ciently lower than the Coulomb exchange energy but still 
finite, the U(l) phase fluctuations are described by the 
following classical action || ^| 



'XY 



J d 2 x (V<5 



(6) 



The unbinding of the vortex-antivortex pairs triggers the 
KT transition from the pseudo spin XY quasi-long-range 
ordered phase to the disorder phase at B c . 

Vortices and anti-vortices are charged || and they 
can contribute to electrical transport. The reason why 
they carry a charge can be understood on the honeycomb 
lattice as follows [Fig. [|(a)]: A defect in the Kckule or- 
der can be visualized as a A or B sublattice site that 
is not dimerized with neighbors, and hence supports a 
midgap state (zero mode). The presence (absence) of 
an electron on such site makes the Kekule vortex or 
antivortex positively (negatively) charged. The pseudo 
spin is pointing T z = +1 (—1) at the vortex core while 
T z = —1 (+1) at the antivortex core. Hence, the charge 
and currents generated by vortex excitations are given 
by 3,i = (T z /ir)e^x9 v dx<p, where /i = 0,x, y. 

In the XY ordered phase where the vortex and an- 
tivortex are bound, it is energetically favorable for the 
T z of the pair to point in the same direction. As a re- 
sult, the bound vortex pairs are charge neutral. On the 
other hand, a charged vortex-antivortex pair can be in- 
duced and pinned by charged impurities. As increasing 
the number of vortex-antivortex pairs, by increasing the 
impurity density or by decreasing magnetic fields, the 
binding interactions between vortex and antivortex are 
screened, causing the KT transition. 




FIG. 2: A typical local XY pseudospin configuration (T x ,T y ) 
represented by arrows. The z-component is represented by a 
grey plot with black as +1 and white as —1. The number of 
orbitals per valley and per spin is — 50, and the number 



of impurities is Mn 



20. 
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In the KT disordered phase (B* < B < B c ), the vor- 
tices are unbound, and their diffusion gives rise to a con- 
ductivity given by a oc n vtx p v t x where n vtx is the den- 
sity of vortices and /i vtx the mobility. This vortex con- 
ducting mechanism is a two-dimensional analog of the 
soliton conduction mediated by charged defects (domain 
walls) in polyacetylene |3C| . In the KT disordered phase 
"vtx ~ l/£ 2 where £ is the KT correlation length^). 
Since £ oc e a ^ Bc ~ B , this gives rise to the KT-divergent 
resistivity. This argument closely follows the one used by 
Halperin and Nelson in analyzing the behavior of the elec- 
trical conductivity of a thin film superconductor above its 
KT transition |5l]]. Indeed, our situation is dual to theirs. 
In Rcf. j3lj the Cooper pair (charge) current exerts the 
magnus force on the vortices and, through the finite vor- 
tex mobility, induces a vortex current perpendicular to 
it. Since vortex current causes an transverse electric field 
(hence a voltage drop) through the Josephson relation, 
this gives rise to a finite electrical resistivity. In our case 
the vortex is charged, and it is the external electric field 
that induced the vortex (charge) current. Thus our elec- 
tric field plays the role of charge (Cooper pair) current 
in Ref. while our charge (vortex) current plays the 
role of electric field in Ref. As the result, electrical 
conductivity in Ref. |3l[] should be translated into elec- 
tric resistivity p in our case; the finding of a ~ £ 2 in Ref. 
|3lf implies p ~ £ 2 in our situation. J32| 

The spontaneous inter-valley coherence discussed 
above is very similar to the inter-layer coherence in the 
double-layer v = 1 QHE || [||. However, there are sev- 
eral important differences, (i) The parameter d/is (d is 
the interlayer separation) in the double-layer system is 
replaced by cl/£b where a is the lattice spacing. For the 
current system a/Is <C 1, a regime which has not been 
achieved in the double layer system, (ii) The inter-valley 
coherent state we propose is spin singlet rather than spin 
polarized. (From this point of view, the v = 1 bilayer QH 
system is similar to the v = ±1 QHE in graphene rather 
than v = 0.) Although the (pseudospin) supercurrent 
cannot be directly measured in the inter-valley coherent 
state in graphene, these two facts have advantage over 
double-layer QH systems to observe the KT physics. 

We stress that our proposal is motivated by the ap- 
parent KT-divergent resistance and the highly resistive 
metallic state at v = observed in Ref. 0]. The very 
basis of our proposal should be subjected to further ex- 
perimental scrutiny, by changing temperature, magnetic 
field, doping, and mobility. 
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